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Abstract 

The second order perturbation calculations for gravity wave and Einstein equation for space 
time and matter are presented for the FRW metric cosmological model. While exact equations 
are found, suitable approximations are made to obtain definite results. In the gravity wave case 
the small wavelength case allows nearly locally flat background for obtaining a flt to the WMAP 
data. In the density and curvature case the FRW background is retained for the length scale 
of WMAP. Clustering and inhomogeneity are understood. The gravity wave ripples from Big 
Bang couple nonlinearly and redistribute the modes to higher values of '1' giving consistency 
with the WMAP results. The order by order consistency of Einstein equations relate the second 
order perturbations in the curvature and density and the wrinkles in spacetime caused by the 
gravity wave modes reorganize these distributions. The radiation data of WMAP gives the 
picture of a FRW spacetime deformed and wrinkled consistent with matter distribution to one 
hundred thousandths parts variation. 

1 INTRODUCTION 

The WMAP data at the observable horizon of 300000 years after Big Bang gives confirmation 
of anisotropy and inhomogeneity in the universe. A bubble model during inflation is made with 
composite equation of state. An explanation of the anisotropy with WMAP C{l)vs"r is made 
with gravity wave mixing causing 'wrinkled' space time from big bang to WMAP epoch. A density 
perturbation model within a Newtonian and Einsteinian analysis is presented. The clustering 
model attempts to explain the two fluid phase with cluster formation at different length scales. 
The recent WMAP picture shows that ^and^^ of the order of 10~^.Then it evolves to clustering 
and clumping of structures. There are number of models for this clustering. In this project as a 
continuation of Part A and Part B.Thc density perturbations are evolved till 1 billion years from 
Big Bang when galaxies are formed. Therefore we simulate the process of gravitational clustering as 
process of diffusion running backwards in time i.e, with an anti-diffusion model. A simple model of 
a uniform Gaussian distribution which can collapse to localized sharp distributions simulates the 
process of variation in density. We have constructed a 1-d Gaussian and a 2-d Gaussian depicting 
variation of density on a line and in a plane respectively. Based on this a 3-d Gaussian representing 
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variation of density in volume may be understood.We are actually using Gaussian distributions as 
solutions of a diffusion equation with the following parameters the width of the Gaussian, the center 
of Gaussian and height of Gaussian. With these parameters a well spread Gaussian becomes narrower 
with time running backward and becomes wider with time running forward. Therefore we simulate 
the process of gravitational clustering as process of diffusion running backwards in time i.e, with 
an anti-diffusion model. This simulation is initially applied at the scale (4000 mega parsec) of the 
Universe leading to formation osuperclusters.Then within a supercluster we understand the process 
similarly to explain the formation of clusters. Finally inside a cluster we associate the formation of 
galaxies with a similar argument. 

The recent WMAP picture shows us that during its initial phase there was no differentiation 
between various forms of matter namely dark matter, dark energynormal matter and radiation. But 
as the Universe evolved different forms of matter started differentiating among themselves and 
occupied different locations. To be more specific the redistribution of normal matter was more pro- 
nounced as large lumps of matter started accumulating together to form centers of galaxies which 
are surrounded by dark matter halos. Each halo is five times more massive than the normal matter 
at the center. Note that the normal matter at the center together with the dark matter halo forms 
the respective galaxy .Some of these galaxies came together to form clusters of galaxies and many 
such clusters approached each other to form super cluster of galaxies. The dark energy as well as 
dark matter occupied in the empty space which are called voids. The redistribution of dark energy 
is not as much as that of normal matter. Therefore an early Universe which had a uniform distri- 
bution of various forms of matter ended up with a one in which the various forms of matter got 
segregated at different places. The title of our report precisely calls for an answer to this observation. 

The overall properties of the Universe roughly speaking are close to being homogenous and yet 
telescopes reveal a wealth of detail on scales varying from single galaxies to large scale structures of 
size exceeding 100 mega parsccs.Thc existence of these cosmological structures tells us something 
important about the initial conditions of the Big Bang, and about physical processes that have 
operated subsequently.If the initial Big Bang would have been homogenous then the presently 
observed distribution of galaxies cannot be accounted for considering the age of the Universe. If the 
elementary particles in the Universe started out with a uniform distribution, then purely statistical 
fluctuations would occur on all scales, from which matter would eventually condense. However the 
highly improbable fluctuations of galactic dimensions is extremely slow and takes far more time 
compared to the age of the Universe. Therefore studying models with inhomogeneities is taken 
up. The aim of studying inhomogeneities is to understand the processes that caused the Universe 
to depart from uniform density The second order perturbation calculations for gravity wave and 
Einstein equation for space time and matter are presented for the FRW metric cosmological model. 
While exact equations are found, suitable approximations are made to obtain definite results. In 
the gravity wave case the small wavelength case allows nearly locally flat background for obtaining 
a fit to the WMAP data. In the density and curvature case the FRW background is retained for 
the length scale of WMAP. Clustering and inhomogeneity are understood. 
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2 Bubble model 



Introduction 

A model of the expansion of the Universe by the expansion of a bubble in some surrounding medium 
when the contents of the bubble are heated is made. The equations of evolution of the bubble radius 
thus obtained will be considered as the equations for evolution of the 'scale factor' of the Universe. 

It is hoped that the dynamics obtained helps to shed light on modeling of inflation, dark mat- 
ter/energy etc. 

We also attempt to find the equation of state of the matter in the Universe if it is assumed to 
evolve according to the dynamics obtained. These equations of state are analysed to find hints of 
dark matter, dark energy and the like. 



2.1 General heating of the bubble 

Consider a spherical bubble containing an ideal gas in a surrounding medium of some liquid - e.g. 

water - at a 

pressure and temperature Pq and Tq . Let the surface tension at the interface be S and for the 
gas in the bubble let ^ = a. At an instant of time say t = heat from a power source - like a 
laser beam or a spark plug - is incident inside the bubble causing it to expand. Let us consider the 
dynamics of such a bubble- 
Now for any radius R of the bubble we have for the gas inside. 




(1) 



Prom the conservation of energy for the bubble we have- 



dQ = nCydT + SdA + PdV 



(2) 



We also have for the bubble- 



A = 4TrR^ dA = STTRdR 



(3) 
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For the ideal gas in the bubble we have- 



T=^ ^ aT=E^^^L^ (4) 
uRn uRn 



Using eqn (g]) in (g]) 



dQ = ^{PdV + VdP) + SdA + PdV 
Rn 

= {PdV)a + {VdP){a-l)+SdA (5) 



Now we substitute eqns dD) and ([3]) in 



2 Q 4 2 S* 

dQ ^ a{Po + —){A-KR^dR) + (a - l)(-7ri?^)( -dR) + S{8TrRdR) 

R 3 R-^ 



Dividing throughout by dt and using W{t) = ^ 



/IS 4 29 

W{t) = R ^a{Po + —){inR') + {a - l){-nR'){- — ) + S{87tR) 

i.e - 



R= , ^(j., - (6) 
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2.2 Casel - Constant power source 

We consider a special case of eqn (jU when the power supphed is a constant. So we have - 

W{t) = Wo (7) 

Then we see that we get the fohowing equation which can be readily integrated- 
J AttR [olPqR + + dR = Wo J dt 

to give- 

^{aPoR + 2Sia + 2)) = ^t + K (8) 
3 47r 
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2.3 Case2 - Bubble conducting heat 

Now lets consider a special case, where the bubble conducts heat across its surface to the surround- 
ing reservoir. Assume that the reservoir is large enough so that it is unaffected by the heat absorbed 
or released by the bubble. 

Let the rate of heat conduction be proportional to the area and the temperature difference 
across the bubble surface. So we have - 



W ^Wc ^ -PR^{T -To) (9) 



We substitute for the T's using the eqns ^ ^ and (g]) 



Air 

{PqR^ + 2SR^) (10) 



N 



Using eqns ^ and (fTO|) we get 



which reduces to 



A-jr 

onJrijsf 



Wc = -^^PR^ {PoiR' ~ Ro) + 2S{R' - rD) 



A-jr 

= ~——l3R^R - Ro) {PoR^ + {PoRo + 2S){R + Ro)) (11) 



Using ini) in 



R ^ --l-RjR - Ro) + f t (12) 

3ni?jv {aPo)R+^{a + 2) ^ ' 
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We note that the sudden increase in R{t) is very similar to the inflation scenario in conventional 
Cosmology. 
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2.4 Other models 

We shall consider variants of previously developed models and compare them to existing cosmolog- 
ical models. 

2.4.1 Hybrid model with inflation 



We consider a model where the bubble first expands rapidly according to eqn but then eqn ^ 
takes over. This involves a discontinuous change in the function W{t) from that given by eqn ([9]) 
to eqn ([7|). This can be considered to be analogous to a phase transition that is presumed to have 
happened during inflation. 



2.4.2 Single model with inflation 



We shall use a variant of eqn where Rq is no longer considered a constant but a function of 
time. We note that in such a case R{t) initially undergoes a sudden inflation and then follows 
exactly the functional form given by Ro{t). 

Thus we use the Friedmann evolution as the form for Ro{t) and so we have an initial inflation 
followed by a Friedmann evolution!!! 



3 Gravitational Waves 



Introduction 

Here we shall attempt to solve the equations of General Relativity for the case of small perturba- 
tions in the background metric and obtain wave solutions. 

Calculations are done upto and 2"'^ order in the perturbations. 

We hypothesize that a primordial gravitational wave is responsible for the anisotropy in the 
CMBR as observed by WMAF. Starting from a isotropic primordial gravitational wave, we hope to 
make use of the non-linear coupling in a gravitational wave solution — upto 2"'* order — to evolve 
the wave into the anisotropics seen in the WMAP data. 
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3.1 Equations of General Relativity 

We list here the equations of General Relativity that shall be useful in this calculation. 
Let us consider a metric gik on spacetime; then we have the following - 



• The Christoffel symbols are 

r*/c( = -j^g^'^{9mk,l + 9ml,k — 9lk,m) 

• The Riemann tensor is 

klm — km,l kL,m ~r nl^ km nm^ kl 

• The Riemann tensor has following symmetries 

^iklm ~ ^kilm ~ ^ikml ~ ^Imik 

^iklm ~i~ ^ilmk ~t~ ^imkl ~ 

_i_ 7?^ -1- n 

ikl;m "r ilm;k ~r imk;l — 

• The Ricci tensor is given by 

Rik = R\lk 

• The scalar curvature is 



Given the energy-momentum tensor Tjfc Einstein equation can be written 

„ 1 „ Snk 

^ik — 29ik^ — —^Tik 



Snk 1 
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3.2 Perturbations 



To analyze the phenomena of gravitational waves, a perturbation hik in the metric of spacetime gik 

ed as - 

gik = + hik (19) 



over its "background" value g'"^ is considered as 



It is chosen by convention to "raise" and "lower" indices using the "background" metric g^^ . 
That is - /i^ = g'^^^^^Kk, where - 

miC) (0) rm /onA 



The perturbations in the contravariant metric is evaluated by solving the equation - 

= C (21) 



Let - (7™* = g""'-^-' + f "'-\ Then substitution into above equation gives 



I.e. 



+ /r + r'h^k = o (22) 

Now, upto first order in h we have - 



mk (■'-) urnk 



f 

Substituting this back in Eq.([?D|) we find the next order as 
Similarly we can iteratively find the higher orders of / as - 
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3.3 Ricci Tensor 

The Riemann tensor is given as 



R^klm — {9nm,kl + gkl,nm — 9nl,km — gkm,nl) (24) 



The Ricci tensor is thus 



Rkm — R klm — {9nm,kl ~^ 9kl,nm 9nl,km 9km,nl) (25) 



The Ricci tensor is now computed upto 2"" order in the perturbation h. 



1 fgln^^^ _ ^In _|_ f 9nm,kl^^^ + 9kl,nrn'^^ — 9nl,km^^^ — 9km,nl^^^ 

nm,kl ~l~ hjzl^nm ^nl^km ^km^nl 



Now separating the Ricci tensor into parts of different orders we have - 

Rkm^ ^ — 2^ {^nm^kl" ^ ~^ 9kl^nrrS ^ 9nl,km^ ^ 9km,nl^ ''^ (^^) 

Rkm^ ^ ~ "^9 {J^nm^kl ~\~ ^kl^nm ^ni,A;m ^km^nl) 

- ^/'"' {gnmM^"^ + 5fc;,nm^°^ - 5ni,fem^°^ ' Qkm^nl^"^) (27) 

nm,kl ~l~ ^kl,nm ^nl,km ^km^nl) 
+ Ih'^K (gnmM^"^ + dkl^nJ"^ ' 9nl,kJ''^ ' 9km J"^) (28) 
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3.4 1^* order calculation with a Minkowskian background 



Now we shall now proceed with the calculations with following assumptions - 

1. The unperturbed background spacetime is Minkowskian. 

=^ gik^^^ — Vik = Minkowskian metric 

2. Only terms upto 1"* order in hik or its derivatives are significant. All higher order terms are 
negligible. 



For a Minkowskian background we also have R^kim = 0; so we need to consider only the 1 
order term of the Riemann tensor 
. From Eq. (|27p. the Ricci tensor of 1** order is - 



RkrrS ^ — 2^ ij^nm^kl ~\~ hkl,nm hnl,km ^kni^nl^ (^'^) 



which reduces to - 

1 , 1 



Rkm = -^Ohkm + -j^ihrn^k,! ^ h,k,m + hkl' ,m) (30) 



where □ denotes the d'Alembertian operator: 



□ _ = A — 

dx'-dx"^ c^dt"^ 



To simplify the equation, we make use of the gauge freedom as observed in Eq.s (??) and (??) 
We impose on hik the following additional condition- 

i'lk = (31) 
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This is known as the harmonic gauge or the de Bonder gauge. 

We again note that the above condtions still do not imply a preferred frame of reference. Sub- 
stituting from Eq. (??) we can see that condition (l3T|) allows for co-ordinate transformations as in 
Eq. (??) provided- 
DC' = (32) 



The Eq. ^ gives us- 



Using this we get- 



And Eq. (??) reduces to 



Rkm = ^n/lfem (33) 



For the case of vacuum we have 



-Rfem — 



and so Eq. ((33|) gives us the usual wave equation- 

nhkm = (34) 
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3.5 2"*^ order calculation 

The wave equation obtained in the previous section is linear. As a resuh there is no self-interaction 
in the wave and so one cannot hope to evolve an isotropic wave into an anisotropic one. With this 
purpose in mind we here attempt to calculate the 2"'^ order corrections to the above equations. 



The second order Ricci tensor is given by Eq. 

RkvrS ^ 2^ {hnva^kl ^kl^mn ^nl^km ^/cm,n/) (^^) 

Thus the wave equation in vacuum reduces to - 

'^/cm — ^ {^nm.kl ~^ hkl^rim ^nl.km ^/cm,n/) (*^^) 

4 Perturbations 



To analyse the phenomena of gravitational waves, a perturbation hik in the metric of spacetime gik 

ered as - 

gik = .gifc*-°^ + hik (37) 



over its "background" value g^^^ is considered as 



It is chosen by convention to "raise" and "lower" indices using the "background" metric g^^^^ . 
That is - h'^ = where - 

^™(o)^^^(o) (38) 



The perturbations in the contravariant metric is evaluated by solving the equation - 

g'^'g^k = (39) 



Let - = g"^'-^^^ + J""**. Then substitution into above equation gives 



K + fk + nhk = (40) 

Now, upto first order in h we have - 



/ 



mk f-^) unik 
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Substituting this back in Eauation (|40p we find the next order as 



Similarly we can iteratively find the higher orders of / as 



ymfc(^) ymz(^ ' jmk^^) hJ^^ (4f) 
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5 Christoffel Symbols 

The Christoffel symbols in terms of the metric tensor are - 

^km ~ 2^ {Qmifk "I" 9ki,m 9km,i) 

where the 'comma' (,) denotes partial derivatives. 
Now, upto 2"'^ order in the metric perturbations - 



So the Christoffel symbols of different orders are - 

r" (0) ^ 1 m(0) („. . (0) , (0) _ M\ 

km — 2 \yim,k ^ y%k,m ykm,i J 

(2)_l,j,p„/ ,(0), (O)_0i. -(0)^ 

km — 2 P yyim,k ^ yik,m, ykm,,i J 

2^ (^jm,fc "I" f^ik,m f^km,i) 
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6 Harmonic Gauge 



To find specific solutions of the Einstein equation, the freedom in its general covariancc is restricted 
by imposing certain additional conditions on the metric called gauge conditions. 

The harmonic gauge, generally used in analysis of gravitational waves, is specified by the con- 
ditions - 

= g'^^Tl^ = (46) 



We now write these conditions for different orders of ft, - 

fTn(0) _ fcm(0)pn (0) 



1 fcm(O) »«(0) 

2^ y 



{gin.,k^'^+9ikj'^-gkm/'^) (47) 



I7n(l) _ „fem(0)-pn (1) , ^femCl)™ (0) 



+ {him,k + hik,m - hkm,i) (48) 



lTn{2) _ km('^)-rn (2) , kmi'^)-rn (0) , fcm(l)T^n (1) 
—9 '-km + 9 '-km + 9 '-km 
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7 Ricci Tensor 

The Riemann tensor is given as - 



R 



Mm — 2 [gkLrni + 9mi,kl 9il,km QkmAl) Q-jp {^''kl^^rni ^^km^^l^ 



The Ricci tensor is thus 



Rkm — R^klm — -^g'" {gmn,kl + gkl,mn — gnLkm — gkm.nl) + g'"gjp (^kl^mn ~ ^km^ 



Now in the case of a Minkowskian background, we have gik^^^ = rjik- So Rkm'"^^ 
first order we have - 



Rkm^ ^ — ^ {hkl^rnn ^mn^kl ^nl^kra hkm,nl^ 



— 2 {^k,nm + ^m,lk h^km + ^hkr. 



Also from Equation (|48p the harmonic gauge condition becomes 

hk,i - \h,k = 



and thus we have 



which in the first order approximation reduces to the ordinary wave equation. 



Also the second order Ricci tensor is 

1 



h {hkl^mn H~ ^mn,/c/ hkrUjUl hnl^km) 



^k,l ~^ ^l,k ^kl'''^ {h^i^n ~^ ^n,m hmn^ 

{K,m + ^Lk - hkm^) (hl^ + - Kr 
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Using the harmonic gauge as in Equation. (|53p we get - 
and the second term reduces to - 

Now we use the condition of Ricci-flatness Rkm = and thus we get the nonhnear wave equation 

'^fcm — h (^/c/.mn ^~ ^7nn,A:/ ^km.nl hnl.km) 

— 2 (^P*:'" + ^p,k ~ ^k,p) (^m,ri + ^n,m ~ hrnn'^) (56) 
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8 The Wave Equations 



Now lets us assume that the metric perturbation h takes a form similar to the linear solution of a 
transverse wave traveling along the Z - axis. So we can write hik as - 



hik = 



f k 

a 

b 

\ 



\ 

b 

-a 

y 



(57) 



where a, 6, k are function of x'^. 



For the case of Minkowski background, the zeroth order harmonic condition is satisfied identi- 
cally. We thus, evaluate the harmonic gauge conditions of 1** order - 



0,n 



2,n 



Thus for difi^erent values n = 0, 1, 2, 3, we get - 

2 dt 

dx dy 2 dx 

^2(1) 9b da Idk ^ 
dx dy 2 dy 

„3(i) 1 dk 
^ =29^=° 
These equations will help simplify the wave equation further. 



(58) 
(59) 
(60) 
(61) 



The wave equation can now be written for each component of hik- For fi = 1,2,3, we have 
ho/j, = 0, and using the first and last harmonic conditions we have - 

= h^^ (/l^l.Ol — hiiflu) + h^^ (/i^2,01 — hi2flfj) + h^^ (/l^l,02 — /l21,0/:i) + h^'^ (/l/i2,02 — h22fin) 

/ Kl (- v°) + Kl (- V°) - K,i Ko) - K,2 {hlo) \ 
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Now for = 1, we get - 



f da db\ d f da db\ 

dt \dy dx J dt \ dx dy J 

/ dk da dk db\ / da da db db\ 

\ dx dt dy dt J \ dx dt dx dt J 



On simplification using the harmonic conditions, we have 



da ( „db da\ db f^da db\ 

mydi^d^J^diyd^^d^J ^ ' 



Similarly for = 2 we have- 



db f^da db\ da f^db da\ 

di\d^^d^)^di\d^~d^) ^ ' 



and for = 3 - 

32 „ q2 



d'^a d% 1 f dada db db 
^dtdz ^ dtdz ~ ~2 ['dtdz didz 



(64) 



For fi = 1,2, 3, hsfj, = so we have - 

= h^^ (/l^i,3i - /lii,3^) + (/i^2,31 - /121,3m) + (^m1,32 " /112,3m) + ^^^^ (V,32 " /l22,3,^) 



2 - 



/^la (/^i.i + hi, - h^^,) + h\, (hi, + h 
hh {K,2 + - K2) + his {Ka + - h'^2) 



\ 



Now for /Lt = 1 
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Using the harmonic conditions we have - 



Similarly for = 2 



and 12 



da da 
dz dx 

da da 
dz dy 



d^a , d% 



dz^ 



dbdb 
dz dx 

db dh 
dz dy 



da 
dz 



db' 
dz , 



Similarly writing the components explicitly wo got for hn 
□a = -k 



\dx^ 



da 
dx 



da 
dy 



2 ^ 


d^a 


d^a 


dxdy 


dx'^ 


dy^ 















da 



dk 
dx 



db' 
dz 



For hi2 = b - 

nb = -k(-^ —\-b(2^^ 

\ dxdy dt"^ J \ dxdy dx"^ dy"^ J 

\dxj \dyj \dxj \dx J \dy J \dy J 2 \dx J 



For /i2i =b we get the same equation as above as is to be expected. 
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For /i22 = —d - 



□a ■■ 



^ d% d'^a d'^a 
dxdy dx^ dy"^ 



+ 



da 
dx 



+ 



da 
dy 



db 
dy 



dk 
dx 



+ 



db 



daV 
dt) 



dk 
dy 



(70) 



Finally for /iqo = k we get 



ak = 



k 



(d 



9^ 
dy'^ 



2b 



d^k 
dxdy 



+ 



dk 
dx 



dkV 

dy) 



- 2 



d'^a , d% 



+ 



(71) 



9 Solution 

We can find a solution by assuming that k = constant and a = Au{z, t) and 6 = Bu{z, t) where 
A,B = constants. For this case the most of the equations identically vanish and we are left with - 

, d'^u ( du\ ^ ( du 

°^ = -^9i^ ' Kdt) =U 

d'^u 1 / du\^ d'^u 1 / du\^ d'^u 1 / du\ ( du^ 

u—TTT = { ^ \ ; U-— = -- — ; u- 



dt^ 2\dtJ ' dz^ 2\dzJ ' dtdz 2\dtJ\dz 



The second equation from these has two independent solutions - 

u{z,t) = Uj^{z + t) ; u{z,t) = U-{z — t) 
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Using these the LHS of the first equation is identically zero giving k = and the other three 
equations all reduce to - 

uu" + ^u'^ = (72) 



where the prime " ' " denotes differentiation hy X± = z±t for u± respectively. 

To solve this we write u' = v and thus u" = and get - 

dv 1 2 
uv— + -V =0 
du 2 



Assuming the v ^ - 

which is readily integrated to give - 

where C is the integration constant. 



dv 1 
u— = —-V 
du 2 



V = Cu-'/^ 



Further using v = we get the solution as — 

u±{z,t) =C{z±t + 6f^^ (73) 
This can be written in the spherical polar coordinates as - 

u±{t, r, e,(t>)=C {rcose ± t + 5f'^ (74) 

As any other function this can be expanded in terms of the spherical harmonics to give the 
various multipole moments for this wave. 



u±{t, r,0,cp)=J2Yl ^'J' (^^) ^" ' 



OO 

EE 

;=o m=-l 

2tv 



where the ji{x) represents the spherical Bessel functions of order I and A; = is the wave 
number of the gravitational perturbation. 

This analysis assumes the existence of a flat background geometry. This can be assumed to 
be valid in a FRW expanding universe provided that the wavelength of the gravitational waves is 
small enough (equivalently large k). Thus the above analysis represents the large k limit of the 
full-calculation. 

To do the full calculation in the FRW background a computational code was written in Math- 
ematica and the Ricci tensor computed for a second order perturbation in FRW space-time. 
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10 Fitting WMAP data 

To model the anisotropy in the CMBR we can expand the fluctuations in terms of the orthogonal 
spherical harmonic functions as - 

oo +1 

^ ^ CiYriBA) 

1=0 m=-l 

where the Y^"^ can be written in terms of the Associated Legendre Polynomials as - 



vr {OA) = (_i)^^ 2; + ia^H)! ,^^^ ^^^^^^ 

To get the coefficients C; we can sample the WMAP plot at various points of interest and 

generate a fit to the curve. 
One such fit produced is - 

l{l + l)|a|2 = l3000exp(-^) J4(^) + 2500 

where J4,{x) represents the Bessel function of order 4. 
Thus the form of C; obtained is - 



Ci = 



13000exp(-3^) J4(^) + 2500 
1(1 + 1) 



Using this form for C; the anisotropy can be plotted as 3D polar plot to get a picture of the 
various intricate patterns in the gravitational wave spectrum. 

Such a plot was made for a sample of 8 points taken from the WMAP plot. 

This represents only a coarse-grain model. For an accurate model more number of data points 
would have to be sampled resulting in a better fit to actual data. 

Though a linear analysis would have given a ripple spectrum but there is an explicit mechanism 
to generate higher-1 modes in a second order calculation. 

The nonlinear coupling, through product of first derivatives, in the gravitational perturbation 
can result in generation of higher harmonic modes in the perturbation. Beginning with lower modes 
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/ = 0, 1,2 one can generate the observed anisotropy caused due to the nonlinear couphng of the 
wave as seen in the nonhnear equation. 

This nonhnear couphng can transfer energy from lower l-modes to higher 1-modes. This results 
in an anisotropic wrinkled space-time. This can be linked to the density perturbations as observed 
in CMBR and WMAP as done below. 



11 Kinematics of distribution of matter through a simple 
model 

11.1 Diffusion equation model 

In our context the word diffusion applies to motion of particles of matter matter may be normal 
matter, dark matter or dark energy. Let us try to solve the diffusion problem in one dimensions, 
the method for higher dimensions follows similarlyDivide the space into thin parts. The process of 
diffusion can be considered effectively as the problem of random walk which deals with the statistical 
progress of a particle equally likely to move forwards or backwards. At any instant of time t half 
the particles in a part at the position x diffuse to the part at x + Ax and half the particles diffuse 
to part &i X — Ax. The variables x and t being continuous require a population density function 
n(x,t);its meaning is that n(x,t)Axis the population occupying a short section Ax. Hence we get 

1 1 

-n(x + Ax, t) + -jn{x — Ax, t) = n(x, t + At) 

Subtracting n(x,t) from both sides, it is supposed that the difference is very small compared to 
n(x,t) itself, we get 

=> ^ [?^(x + Ax, t) — n(x, t)] — ^ [n(x, i) — n(x — Ax, t)] = n{x,t + At) — n(x, t) 




dn 

'at 



At 



-(Ax)2— = l — \At 
2^ ^9x2 \dt/ 

The final step is to assume that, with a large number of particles, Ax and At can be reduced 
indefinitely (which makes the above approximations exact) but keeping the ratio D = ^^^^^ 
constant. Thus we have the diffusion equation as 

S^n / dn 

which means that the random diffusing process which we envisage does indeed lead to the diffusion 
equation. A special solution of diffusion equation is the Gaussian which is 



n(x, t) = {nDt) 2 e 



27 



The above equation gives the distribution of number of particles with respect to position at a given 
time t. If there was very distribution at xq at time to then we write the Gaussian as 



As mentioned in the introduction we are trying to model the redistribution of matter in the Universe 
as a process of diffusion running backwards in time. Looked at the scale of the Universe the various 
points of accumulation of normal matter correspond to super-clusters. At the scale of super- 
clusters points of accumulation of normal matter correspond to clusters and the scale of clusters 
to galaxies. Firstly we explain the clustering and clumping of normal matter at the scale of the 
Universe. Similar arguments are valid at the scale of super-clusters, clusters and galaxies. Now if the 
process of clustering of normal matter is diffusion running backwards in time and equation [75] is a 
solution of diffusion equation and to denotes the present time then xq would be the position of a 
typical super-cluster. But equation [75] would mean that all the matter is concentrated at the point 
Xq at present time. This would mean that a super-cluster has no size which is wrong. Therefore we 
add a — ct^ term to equation 1751 where a is the standard deviation which is an estimate of the size of 
super-cluster. For an Gaussian a corresponds to half of the width of the Gaussian at half maximum 
which is 65% of the total width of the Gaussian. For a Gaussian representing a super-cluster the 
radius of super cluster is roughly 3cr as 3cr = 95% of all matter is within the cluster entity of this 
radius. 

Now having modeled the distribution of normal matter through a diffusion equation we can 
account for the redistribution of dark-matter and dark-energy also in a similar way. The most 
important point to be noted is that the Gaussian for different forms of matter will have different 
values of Diffusion coefficient D and different values for standard deviation a. This is the main 
reason why different forms of matter differentiate among each other and ended up at different 
places in the Universe. 1-d as well as 2-d plots were made using Matlab of density of normal matter 
v/s position for various times at the scale of Universe. From these figures one can see how a uniform 
distribution of normal matter ultimately led to the formation of super-clusters with sharp densities. 
The formation of clusters within a super cluster and of galaxies within a cluster can be explained 
similarly. Note that in the real diffusion process of a gas the series of pictures would be from sharp 
to wide gaussians. The radius of the Universe is AOOOmega — parsecs lO^^m.In the 1-d plot 1 
unit=10 mega-parsecs and in 2-d plot 1 unit=100 mega-parsecs were used. 




(75) 
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We now draw a rough analogy between the distribution of various forms of matter in the Universe 
and condensation of air. The rough mass composition of air in our atmosphere is 70% of N2 27% 
of 02,1% of CO2,0.5% of H2O and 0.1 % of CO. As the temperature keeps dropping different 
gases start to condense. The process begins with the condensation of water-vapor. At various 
points in air one finds formation of droplets of water at different places as we go closer to 0°C. 
Condensation is nothing but bringing water molecules which were far closer to each other until 
sufficient accumulation takes place as the phase changes from gas to liquid. As the temperature 
goes down further to —AO'^C the CO2 in the atmosphere also condenses and similarly various gases 
condense at various stages. Now the redistribution of various forms of matter in the Universe is 
analogous to condensation of air in the sense that various forms of matter accumulate at different 
points in space. Note that the mass composition of our Universe is 69.7475% of Dark Energy present 
in voids, 29. 2475% of Dark matter out which 2.5% is in galaxy halos and the rest is in voids, 0.5% 
of normal matter and .005 % of radiation present everywhere in the Universe. 



Name 


Number 


Size 


Mass 


Universe 


1 


order of 10'''' 
meters (around 
4000 mega par- 
sees 10^ hght 
years 


10^*^ solar 
masses 


Supercluster 


roughly 1000 


400 mega par- 
sees 


10^^ solar 
masses 


Cluster 


256000 


40 mega parsecs 


10^^ solar 
masses 


Galaxy 


10^^ galaxies in 
the Universe 


4 mega parsecs 


lO'"* solar masses 



Now we try to establish a criterion to find out whether or not a given set of cluster of galaxies 
form a super-cluster. Let the size of each grid be 1 unit x 1 unit. Now from the figure one can make 
out various patterns of super-cluster of cluster of galaxies. But can we try to obtain a quantitative 
result in order to define a super-cluster. For this purpose consider any pattern from the figure. For 
example we consider the pattern of clusters with coordinates: 

(28,30);(27,32);(30,32);(26,34);(29,34);(31,36);(29,37);(31,38);(29,40) 

Since different clusters are separated from each other by different distances we need to obtain a 
quantity which represents the average value of all distances. Now if the evaluated mean is less than 
some fixed value then the given group of clusters form a super cluster. If (cci, jji) and {xj,yj) are the 
coordinates of the ith and jth cluster then the following quantities may be evaluated for a set of n 
clusters : 
1. 

2. 

\ "C2 



= 10.203 



5.102 



29 



j=i+i 



n 



"C2 



30 



Sl.no 


positions of cluster 


Ax i.e, X 
separation 


z\v i.e V 
separation 


(Ax)^ + (Ay)^ 


((Aa;)^ + (Aw)^)' 


1 


(28,30)(27,32) 


1 


2 


5 




2 


(28,30)(30,32) 


2 


2 


8 


2^2 


3 


(28,30) (26,34) 


2 


4 


20 


2\/5 


4 


(28,30)(29,34) 


1 


4 


17 


717 


5 


(28,30)(31,36) 


3 


6 


45 


3^/5 


6 


(28,30)(29,37) 


1 


7 


50 


5\/2 


7 


(28,30) (31,38) 


3 


8 


73 


V73 


8 


(28,30)(29,40) 


1 


10 


101 


Vioi 


9 


(27,32)(30,32) 


3 





9 


3 


10 


(27,32)(26,34) 


1 


2 


5 


2\/5 


11 


(27,32)(29,34) 


4 


4 


8 


2\/2 


12 


(27,32)(31,36) 


4 


4 


32 


4V2 


13 


(27,32)(29,37) 


2 


5 


29 


\^ 


14 


(27.;-')2)(;U.38) 


4 


G 


52 


^52 


15 


(27,32)(29,40) 


2 


8 


68 


2\/T7 


16 


(30,32)(26,34) 


4 


2 


20 


2V5 


17 


(30,32)(29,34) 


1 


2 


5 


\/5 


18 


(30,32)(31,36) 


1 


4 


17 


%/l7 


19 


(30,32)(29,37) 


1 


5 


26 


V26 


20 


(30,32)(31,38) 


1 


6 


37 


V37 


21 


(30,32) (29,40) 


1 


8 


65 


2\/65 


22 


(26,34)(29,34) 


3 





9 


3 


23 


f26 34H31 361 


2 


2 


8 


2V2 


24 


(26,34)(29,37) 


3 


3 


18 


3a/2 


25 


rafi 341 C31 38") 


5 


4 


41 


2\/41 


26 


f2fi 34V2q 40'! 


3 


g 


45 


3\/5 


27 


C2q 34'l('31 36") 


2 


2 


8 


2v'2 


28 


^29 341 r29 371 





3 


9 


3 


29 


f29 341('31 381 


2 


4 


20 


2\fh 


30 


(29,34)(29,40) 





6 


36 


6 


31 


(31,36)(29,37) 


2 


1 


5 


x/5 


32 


(31,36)(31,38) 





2 


4 


2 


33 


(31,36)(29,40) 


2 


4 


20 


2V5 


34 


(29,37)(31.38) 


2 


1 


5 




35 


(29,37) (29,40) 





3 


9 


3 


36 


(31,38)(29,40) 


2 


2 


8 


2\/2 
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12 Density perturbation using FRW metric model 

We discuss the model based on the Friedman-Robertson- Walker (FRW) model. The FRW metric 
is given by 

ds^ = c^dt^ - a^{t){dr^ + r'^{def + sin^ Q{d^f) 



Einstein's field equation 



where i?^,^ is the Ricci tensor, is the Einstein tensor, R is the Ricci scalar and T^^ is the 
stress-energy tensor. With the above two equations get at the following two equations H 



— + — - 

a 



J- n 



where k is the curvature constant. 

Assuming that the Universe is filled with Pressure less dust then the stress energy tensor T^^ 
becomes 

^ p ^ 


V / 

Hence assuming a pressureless dust-dominated Universe and taking the curvature constant to be 
in the above two equations we get the following two equations. 



2a a2 
- + - = 
a 



a 



w 



(76) 
(77) 



Solving equation [77] we get 

If to represents the present time then 



Po 



3c4 / a 



We now solve equation [77| by rewriting as 



to 

SvrG 



tt2 "O 



(78) 



^Reference 1 Chapter 4 
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a2 



da — 



Hoa§ dt 



2. 3 s. 
^ a2 = -Hoagt 

Assuming a=0 at t=0 the arbitrary constant that arises out of integration is set to 0. 
Putting 

Ho^i^ (79) 
1 3 2 1 

Putting Hq — ^ where to is the present time. 



3to 



=>a = ao[^y (80) 



-1 

One may as well check that this solution of a satisfies [76] by substituting d = |ao (^t~^ ^ ^ ^^^^ 
As we are considering a model of Universe in which mass (M(t)) is conserved we have 





fdipit)V{t))\ 



= 



\ dt J 

where p(t)a,iidV{t) are the density and volume of the Universe respectively. 







V{t) = bR^ where R is the radius of the Universe which is a function of time and b is a constant. 
Put p{t) = p then 



If^-^and%)^^then 



dt — "•'-^'^ dt 

3pR = -Rp 



Since | = ^ 

R a 



p —3R —3d 
p R a 
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Finally 

P- = ^ (81) 
p a 

p —'SR —3d 
p R a 

Conservation of mass would imply that 
Using 

a _ R 

Go Ro 

Substituting ' function of time from equation [5D] we get 

(82) 

Till now we have discussed equations and solutions which are valid provided our Universe is 
isotropic with respect to space i.e, if the scale factor 'a' were function of time alone. The FRW 
model is conventionally done in spherical-polar co-ordinates. Now if isotropy and homogeneity is 
broken we modify a(t) so that f| 7^ and ^ 7^ O.For a 3 sphere we write down a{t, 9, (f>). 

a{t,e,<j)) = ai{t) +a2{t,e,<j)) (83) 

Similarly for the matter part we write down p{t, 6, cj)) 

p{t,0,cj,)^pi{t)+p2{t,0,(f>) (84) 

If P represents the pressure and T the temperature at WMAP time we have ^ ~ ^ ~ ~ 10^^ 
which increases with subsequent time. Hence ^ ~ ^ ~ 10^^ at WMAP time {10^^ seconds. So in 
the equations I76I77I and [81] we include the effect of 02 and try to obtain the new set of equations 
with ^ < 1 , y < 1- 

The detailed structure formation on the Universe as seen at z (doppler shift) =0 to 1 with 
voids, cluster of galaxies and various kinds of fluctuations requires super computing models to put 
observed data from HST, SDSS etc which is outside the scope of this project. We solve for density 
fluctuations from WMAP i.e,z=1000 till z—5 and density fluctuations from z—5 to z=l. Recalling 
that 

2 



1 + z 



a{t) V t 



Now combining [75] and [77] we get 



-+(^)4-0 (85) 
a V 3c^ / a-^ 
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Introducing inhomogeneity and anisotropy in a = a{t, 9, (p) as explained above we have 

ai + a2 30^(01+02)'^ 
Note that oi = ai{t) and 02 = a2(i, 9, 4>) 



= 



1 + - ) + ( 1 + - 1 ' = 



oi \ fli / 3c'*a? \ ai 



In the binomial expansion (^1 + f^^ and (^1 + f^^ we neglect the terms having the power of 
02 greater than 1 

^ 2(ai +02) A _ a2 \ ^ STrGpoag / _ SoaX ^ ^ 



Approximating '^^'^^ to we get 
aiQi 



2di 2dia2 ^ 2a2 ^ STrGpoOo 3*^2 STrGpoOo 



Using equation [55] we have 

2di STrGpoflo \ "1 _ n 



Therefore 



02 J— -3 02 = (86) 



ai 

Now since ai plays role of a in equation [85] and using equation [78] we get 

oi ^ -Hi /a| 
fli 2 ya^^ 

Using the above expression and equation [75] equation [51] can be written as 

.. Hi fal\ 3Hi fal\ 

a? 



^ (i2 - i/oM ^ «2 = 



1 

Since ai plays the role of 'a' in equation 1801 we know oi as a function of time and using equation 
[79l to substitute fori/o we get 

- ^ = (87) 
Plotting a graph of a2 v/s time gives (02 = 02) 
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Introducing inhomogeneity and anisotropy in p = p{t, 6, (p) as explained above we have 

-3(a'i + a2)(pi + P2) 



Pi + P2 

Note that pi = pi{t) and p2 = P2{t, 0, (j)) 



ai + 02 



. . -3(ai + a2)(/9i + P2) /, , a2 
Pi+ P2 = 1 



-1 

02 



In the binomial expansion (^1 + j we neglect the terms having the power of 02 greater than 1 

^ . , . -3(ai +a2)(pi +P2) A 02 
^ Pi+ P2 = 1 

. ■ -3aip2 f -o.2\ -3a2P2 -3a'2pi f -a2\ , -3a2pi /-a2\ , „ 

Approximating , , ana to U we get 

ai \ ai J ai ai \ ai J ai \ ai 



Pi + P2 



3aipi 3dip2 3dipi f —a2\ 3a'2/3i 



Using EI] 



Therefore 



ai ai ai \ ai j ai 

pi _ -3ai 
Pi ai 



3dip2 3aipi f a2\ M2P1 

P2 = + — (88) 

a\ a\ \a\ J a\ 

We have simulated the process of structure formation in our Universe as a diffusion process 
running backwards in time. This means that the density of normal matter satisfies the diffusion 
equation. Actually the process of structure formation of normal matter is anisotropic in space. This 
means that the p2 term of the density satisfies the diffusion equation i.e, 

^ = £'V2p2 (89) 
where D is the diffusion coefficient. But from equation [88l also satisfies 

3ai/92 , 3ai/3i ( a2\ 3a2/Oi 

P2 = + — 

a\ a\ \a\ ) a\ 

Now this is analogous to a diffusion equation since this can be rearranged as 

. 3ai/92 3a'i/9i ( a2\ 3d2Pi 

P2 + = — 

ai ai \o-i / ai 

When analogy is drawn between this equation and the diffusion equation we see that the L.H.S is 
analogous to diffusion equation and the R.H.S suggests the existence of a source within the system 
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undergoing diffusion. For the time being let us assume the source term to be absent i.e, remove the 
perturbation in a. Now in the absence of the perturbation term the equation becomes 



3aiP2 .„„x 

P2 = (90) 

ax 

Compare the above equation with diffusion equation in one dimension for simpHcity. Therefore from 
equation [Ml we get 

From the above two equations it follows that 

^ JJlEl (91) 

This means that DV^ acting on p2 is equal times p2-Note that ^ is a constant with respect 
to position .In such cases the variation of density with position is a sine function and cosine function 
or a linear combination of both i.e an exponential function. One must also observe that is just 
3 times the value of Hubble's constant at different epochs. In equation [89] if is expanded in 
spherical polar co-ordinate system as 

= + + v| 

where ; i = r, , ip are in general functions of all the three variables r , 9 , (j) 

If anisotropy and inhomogeneity is taken via P2it, 9, (p) and 02 (i, 9, </>) then diffusion equation 1891 
takes the form 



From equation 1901 we have 



From the above two equations we get 



P2 



M1P2 



ai 



3a'ip2 



ai 

Now the solution for the above equation will be of the form 

I 

We consider the sampling of data of CMBR improving in resolution from COBE to WMAP to 
Planck satellites will give better angular correlations in 9 and (p. So that 02 and p2 the anisotropic 
perturbations of space-time geometry and matter density can be given values at neighboring 9 and 
9 zt d9 and the difference equation converted to a differential equation. Hence a acts on p2 in 
equation [89l This feature is repeated at scale of super-cluster, cluster and galaxy. 
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Consider equation [551 and substituting ai,di using equation [501 and substituting for rho from 
equation [82] we arrive at the equation 



P2 



-2/32 , Po f 2a2 



-(— -3a2|(^ 



t flo 

Considering the equation 1501 substituting ai,ai using equation [501 we get 

-2p2 



This would be the case if in equation [ 



P2 = 
202 

t 



t 



Mo 



From equation [79| one can conclude that 

H = 



2 

it 



02 
02 

O 

o 



2 

O'l + 02 

Ol + 02 



(92) 
(93) 

(94) 
(95) 



From the above two equations it follows that the anisotropic part of Bubble's constant scales in 
the same way as Hubble when diffusion without a source takes place. This is exactly the situation 
during WMAP when clustering has just started. From equation [OH 



P2 



-2(02 , Po f '^a 



t 



-(— -3a2l 



Oo 



If 02 is approximated to in the equation [92] which is the present state of Universe (from 1 billion 
years onwards) where in the change in perturbation is made very small and using equation [80] and 
equation [82] in the equation [92] we get 

2j02 _ Pi f 202' 
t ai \ t 



P2 



This is an equation in which R.H.S contains the anisotropic part 02 and roughly models the present 
state of our Universe. 

Considering the equation [Oil substituting ai,di using equation [501 we get 



D 



d^P2{z) 



-2p2 
t 



(96) 



By inspection p2oexp 



is a solution where Oon is a constant 

,(^)'.' 

The argument on Differentiation 

If we assume that the Universe is filled with dust having then the stress energy tensor T^i, 
becomes 



f p 








' \ 





-p 














-p 





V 








-P J 
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Then equation [5T] gets modified as 

p + 6F a 

The equation of state determines P for each form of matter is 

P = wp (98) 

where w =-1 for dark energy ,w = | for dark matter — | for normal matter and ^ for radiation. 

Therefore for each kind of matter we have Pi — WiPi Hence pi + 3Pi — pi{l + w^). Hence each 
kind of matter will distribute differently due to different values of p and P. Therefore in equation 1971 
the value of 

p + 3P = ^(p.(l + i«0n,) 

i 

where rii is the number fraction of each type of matter. 
Consider equation [551 

3dip2 , Sa'iPi f M2P1 
P2 = + — 

In the above equation pi is the cumulative of all forms of matter and 02 is also the cumulative of all 
forms of matter since it is related to the space time structure and not on the kind of matter. Therefore 
the only the matter dependent property is p2 -Therefore we say that matter does not differentiate 
itself at the level of pi But the differentiation occurs at the level of p2 and p2 is a function of 
position this means that different forms of matter occupy different positions. 

P2 is undifferentiated at WMAP stage and ^ is of the order of 10^^. Each component is 
evolving differently by the time we get to galaxy formation epoch. That means equation [88] is now 
satisfied independently for each component. Therefore the R.H.S of equation 1881 is common but p2 
is replaced by (p2i(l + Wi)ni). where i runs from 1 to 4.i=l represents dark energy;i=2 represents 
dark matter;i=3 represents normal matter;i=4 represents radiation. The analogy of equation l88l with 
the diffusion equation [89] with the Gaussian solution require that each component has a diffusion 
coefficient an dcr consistent with the observed data. Hence clustering and differentiation has been 
explained. 



13 Two fluid clustering analog model of the Universe 

The Navier Stokes equation of fluid dynamics 



dv 
dt 



{v.V)v 



'VP + TfJ^V+ C 



grad{div{v)) 



Consider an incompressible fluid i.e, let div{v) = Therefore the above equation becomes 



- + (..V). 



-VP + rjV^v 



Now at the scale of universe we say that dark energy is a fluid which pushes dark matter and 
normal matter together to form super- clusters of galaxies but once the structures are formed the 



^is still to be confirmed through experiments 
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dark energy it pushes them apart and the super-clusters move away from each other. One must 
zoom in and find similar process happening at various scales. Assuming that the dark energy is a 
non- viscous fluid in our model we put t] — 



dv 



-VP 



If {v.W)v << VP The equation of state for dark energy is P = —p. Substituting this in the above 
equation we get 



dv 



Vp 
P 



(99) 



According to the equation of continuity 



Now p(yv) 



| + v(H = o 



|-l-V(p).^=0 



But from [Ml we substitute W{p) in the above equation. Therefore 



dp 
dt 



dv 



dp 
p 



'V.dv 



(We have taken dark energy to be a fluid whose density is invariant with position) Therefore 

Hp) = — + ^ 

where p 

is the density of dark energy,?/ the velocity of dark energy and b is the constant of integration. 
In case of dark matter the fluid equation is obtained as follows. Consider the Navier-Stokes equa- 
tion. The clustering of dark and normal matter together suggests that VP — and the formation of 
various kinds of galaxies like spiral suggests that the viscous force is not negligible. But the fluid is 
assumed to be incompressible i.e, diviv) = 0. Therefore the equation for dark-matter normal matter 
mixture becomes 

dv 

— + {v.\/)v = r^V^t/ 

The flow of any fluid depends through the boundary conditions, on shape and dimensions of the 
body moving through the fluid and on its velocity. Since the shape of the body is supposed given, its 
geometrical properties are determined by one linear dimension, which we denote by l.Let the velocity 
of the mainstream be u.Then any flow is specified by three parameters, j/,u and 1 where nu = ^ is 

2 

the kinematic viscosity. These quantities have the following dimensions: i/ = , I — cm , u = 
I — cm , u — —. It is easy to verify that only one dimensionless quantity can be 



cm 
sec 
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formed from the above three, namely ^.This combination is called the Reynolds number and is 
denoted by R: 

pul ul 



R = 



V 



In case of the dark-matter, normal matter fluid one can assume the flow to have small Reynolds 
number because of high value of viscosity of fluid. Now the term {v.\7)v) is of the order of magnitude 

2i-.The quantity V^v is of the order of magnitude The ratio of the two is just the Reynolds 

number. Hence the term ((w.V)u) may be neglected if the Reynolds number is small, and the 
equation of motion reduces to 

One can observe that this is nothing but the diffusion equation for the velocity vector. The velocity 
vectors which are functions of position and time satisfy the diffusion equation with diffusion constant 



14 Anisotropy and inhomogeneity of Density in universe 

In standard models of Cosmology like Priedmann Robertson Walker ( FRW) models the universe is 
supposed to be homogenous and isotropic. But as we zoom in to smaller scales the inhomogeneity 
and anisotropy in the form of superclusters, clusters and galaxies appear. 



The formation of the Large Scale Structure takes place due to perturbation in the density 
of matter. Our aim is to introduce a small perturbation in the density of the universe. This 

perturbation is a function of 0, (f) and t. 

The dynamic equations for the FRW-Model have the form- 

" = (101) 




The above two equation's are not independent of each other and from the above we obtain the 
relation- 

R f rV , , 

r = -\r) (10^) 

In equation (28) we put 

p = p{t)^bp{d,4>,t) 
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and 



R= R{t)+5R{e,(t),t) 



where 5R and 6p are perturbations of first order in R and p. 
The equation (28) after solving leads to 

In deriving the above equation the assumption is made that R{t) and p{t) follow Eq.. (28). In the 
above Eq.. i?,(t), R{t) and p{t) can be derived from Eq. (28) as functions of time by making the 
substitution p = ^"^y* . Thus by knowing the variation of 5p with time an ordinary differential 
equation governing the variation of 5R with time is obtained. 



In order to obtain the variation of 5p with time we make use of the equation which gives the 
variation of total density with time. 

^--4 ,104) 



The above equation is derived by assuming that the total mass of the universe is constant. 
Substituting ^ from Eq.. (28) in terms of p leads to- 



p=-3p3/y!^ (105) 
In which after substituting p=p{t) + Sp{0, </>) and assuming that p{t) follows Eq.. (31) we get 



Sp 9 SnGp 

Tp = -2^^ 

From Eq.. (33) the time variation oi 5p can be obtained and thus the time variation of 5R can 
be obtained. 

15 Distribution and differentiation of various forms of mat- 
ter in the Universe 



Consider the differential equation of a2v/st 



42 



The solution of the above equation is + C2t^ where ci and C2 are integration constants. One can 
see that this function attains an extremum at — 

ci 

Supposing the initial conditions were (i.e a2 = at t = we get the solution as 



a2 = c{0,(j))t3 
Consider the differential equation p2v/s t 

p2 = + — 

ai ai 

Substituting equation If 071 in the above equation along with 



(107) 



3a2/Oi 
ai 



and 



Therefore 



P2 = 



P 



ai = ao I — 



P2 i^-Aa^ii-Y]^ 



ao 



Pi = 



2p2 ^ 2c{e,c^)potl _ 4c(g,0)pot| 



3^ (c(e,0)it)^# 



ao 



(A)' 



We have — f3.7 x fO^ years , po = 2.f f x fO^^^ kg/m^ and ao = f Hence 

2p2 i.blx 10^»c{9,(t)) 



P2 



P2 = exp 



t ti 
4.5f X 10^^c{e,, 



P2 = t 

P2 
P2 



4.5f X fOi^c(^,(^) I [t'-i^ 
4.5f X fO^^cO 



exp 
i\ dt + K 



jdt] ] dt + K 



t— 



6.765 X fO^*c(6i,0) t 



If we substitute the initial conditions at t=WMAP i.e — 

ai+a2 



Kt-' 
= 10~^ we get 



t3 



+ C2t3 



f + C2tl+ao(^)' 



= 10" 
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2.00276 X W^C2 + 4.72702ci 7.82709 x 10* (108) 
However if we substitute 02 to be -^f + C2t^ in the differential equation of p2 (if one is not sure 

t3 

about the initial conditions of 02) we get 



. _ _2p2 ^ 2cipoto^ ^ cipo^g ^ 2c2Pot§ _ 4c2pot(f 
t flo^^ ao^^ ao^J flot^ 

We have to = 13.7 x 10^ years , po = 2.11 x 10^^^ kg/nr" and aa — 1 Hence 

~2p2 6.76502 X lOi^ci (-4.51001 x IOI8C2) 



P2 



Pi 



t ti 
2p2 _ 6.76502 X lOi^ci (-4.51001 x IOI8C2) 

~" ^ ^1 



Call 2 as A:i,6.76502 x lO^^Ci as k2 and— 4.51001 x 10^^C2 as k^. Therefore the above equation 
becomes 

. fci k2 kz 

^^ + T = F + Ii 

The above equation is a linear first order differential equation whose solution is 

where k is the constant of integration. Substituting the values of fci, k2 and k^ respectively we get 

P2 ^ (-6.76502 X 10^**cit"3) + (-6.765 x lO^^cai^) + kr'^ 

where k is the constant of integration. 

We know that at t = WMAP = 300, 000 years 

= 10-5 (109) 

Pi + P2 

P2 sXt = 300,000 years is calculated to be -7.9724 x lO-^^ci - 33.778c2 + 1.116 x 10"^^ and pi 
at t = 300, 000 years is equal to 4.40029 x lO^^o Therefore at t=WMAP substituting pi and p2 in 
[TUgiwe get 

3.3778 X 10% - 7.9724 x 10"^% + 1.116 x 10"^^ = 4.40029 x 10"^° 
Now ci and C2 are actually functions of 9 and Therefore imagine ci = af{0, (/)) and C2 — Pg{0, (p)) 
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16 Calculation of Einstein tensor for the modified Fried- 
mann Robertson Walker metric 



Let us call a[t] as R[t]. The modified Robertson Walker metric is given by 
{dsf = {dtf-[R{t)+b{t,0,(P)]\drf-[R{t)+b{t,O,(l))]\rf{def-[R{t)+b{t,e 

where b{t, 9, 4>) is the perturbation introduced in R{t) Then the components of the Einstein Tensor 
for the modified Robertson Walker metric which introduces inhomogeneity is 

Below R is a function of t and b is a function oft , 6 , ^.Note that for each of the terms below the 
term in the first bracket is nothing but the terms of the Einstein tensor for the Robertson- Walker 
metric. Also in evaluating the Einstein tensor a binomial expansion was made and terms involving 
higher powers of b or the derivatives of b have been neglected. The complete Einstein tensor was 
evaluated using Mathematica. 



17 Einstein tensor for perturbed FRW metric 



00 



Gil = 



(rR)^ 



G22 = 



{rRy 



r'^R^ 



[3{rRR')^ + b (6R{rR'f - 12{rR' fR 



+^{-2Rcom) + 



G02 = 



Go3 



i?2 

2 

R? 



-R 



-R 



= 

+ R 



^{6{rRfR') 



d-'b 



{-2Rcosec^{0)) 



db 
80 

db\ 



/ d^b Y 



+ R 



G 



10 



[rRR'f + 2r''R^R"j + b (2R{rR'f + 6{rRfR" + ^ (^{rRR'f + 2r''R^R"^^ 
+fA2{rRfR:) + m-Rcom) + ||(2r2i?3) + ^{-R) + ^^-Rcose^0) 



G 



Of 

db 
90 



12 



-1 

rR 



1 

rR 



G20 = 



2 

R^ 



(rRR'f 



2r^R^R 

{2{rR'fR') 



- b {2R{rR!f 



R 



\dtd0) 



db 
80 



\{-Rcot{0)) 



6irR)^R" + ^ [{rRR^f + 2r^R^R"^^ 

dt 



(2r2i?3) 



i-Rcosec^0) 
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^23 — — 



G30 = — 



2 



R 



1 

rR 



Gm — — 



-) (- 

del) ' ^ 



-cote) 



\dtd(t>) 



d0d(t) 



^33 = ;^ 



(rRR'siney + 2r^R^sin'^eR" 
+h{2R{rR sindf + &R"{rRsinef - -| (rRR sineY + 2r'^R^sin^0R" 



+ §{2R{rR' sin9)^) 

If we consider the standard stress-energy tensor (which is dust with pressure) in which the off 
diagonal terms are zero then 5(t, 9,(f>) =0 which removes the inhomogeneity which is wanted in 
the first place. So one must think of an appropriate stress energy tensor which when equated to the 
Einstein tensor gives a solution for b{t, 0, <j)) 

The off diagonal terms in Einstein's equation are set to zero in this approximation as case (1), 
and in the right hand side, the off diagonal elements of the stress energy tensor are set to zero. 
Consider only the diagonal terms. 



^{2r^R^sm^0) 



18 Approximate solution 



Now since ao = 1 



and 



Also 



Also 



R= ai = ao( 

to 







3to 


Kto, 






94 


K.t^ 



b = a2 = c{0, 4>)t3 

gi = r^0A)t-^ 



db^dciMl.i 

80 80 
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Neglecting the higher derivatives of 6 and (j) in Goo we get 



00 



-1 



(3{rRR'f^ + b (-6{rR'fR 



Goo 



{if 



3 r 



\to J 3*0 \to J 



de 

-ix 2 



m-^2Rcom) 



3 



Therefore Goo simpUfies to give 



-1 / t 
To 



{ Scroti 
— (-1 + 



2>tQ \to 



3*0 V *o 



de 



Neglecting the higher derivatives of 6 and 4> in G\\ we get 



Gn = 



{rRf 



{rRR'f + 2r'^R^R" + b (2R{rR'f + Q{rRfR + ^ {(rRRf + 2r^R^Ryj 
+ §{2{rRfR') + §i-Rcom) + ^{2r^R^) 



Gn = 



(*o) 3*0 (to) 



+ 2r^ 



ay 



-2 / ±_ 

9*^ Uo 



2T 3 



9*J V*o/ 



dc(0,4>) 

de 



G 



11 



2\ 2 
3 



28cr2i 3 5c(6', < 



Observation yields 



G 



961 



R 



11 



to 



t \ ^ Scr'^t^ ( t 

cot{e) + —— 



= G22 



Neglect the higher derivatives of 6 = c{0, 3 yields 

Gn = G22 



to 
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Gssdoes not yield a differential equation for b{6, (j)) and hence it is redundant for this approximation. 
Therefore 



r - 

Lrnn — I — 

10 



Gil = G22 = 



4^2 / 






to 









Scroti 



3t 



de 



t3 



to 



ae 



Substitute for t the WMAP time which is 3 x 10^ x 365 x 24 x 60 x 60 and for to 13.6 x 10^ x 365 x 
24 X 60 X 60 seconds. These in the units c = 1 will be measured in meters i.e 

twMAP = 3 X 10^ X 365 X 24 X 60 X 60 X 3 X 10^ meters = 2.84 x 10^^ meters 

to = 13.6 X 10^ X 365 X 24 X 60 X 60 X 3 X 10^ meters = 1.29 x 10^^ meters 
The first term in Goo is far less in magnitude than the other terms. Hence 

Scr'^ti 



-1 / t 



Goo = I — 



Gil = G22 = 



to 

1 / t 

Tn 



'itl 



^cr'^t^ A 
2 cot{e)V' 



t \ ^ dc{e, (j)) 

to 



de 



Now in c = 1 units Gqo = 87rG(pi + ^2)- This includes the first and second order terms. The Goo 
which we obtained is second order, and approximating p2 ~ 0. 
8nGp2 has a very small value compared to terms on the L.H.S. 

19 case 1 



-1 / 1 



r' \to 



— 2 — 2cot{e)t^ ( — 



2,t 







to 



de 



Scroti ^ ^ f dc(e,<P) 

— ^ — 2cot{e)t3 ( — ' ^ 



3ti 



to 



80 



4r2 {twMAp) ' (to) ^ 



tan{e) de 



c{6A) 

d[c{e,ct>)] 

c{e,<p) 



1.024 X lQ-'''^rnn{sec{e)) + ln{g{(j))) = ln{c{0,(t))) 
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where g{(l)) is the constant of integration. 

Now in c = 1 units Gn = G22 = 8nG{Pi + P2). Now the Gn and G22 which we obtained is called 
second order and approximating i.e P2 » then SirGP^ has a very small value compared to terms 
on the L.H.S. Hence we get 



-1 { I 



4crH^ 4 / i \ " dc(9, (j)) 







(twMAp)^ ito)^ [tanie) 80] = (HO) 

d[cie,<i>)] 



=^ 4r" (twMAp) ' (to) ' J tan{9) 09 -- j 

3 X 1.024 X IQ-^^ r"^ ln{sec{e)) + ln{h{(l))) = ln{c(9,(j))) 
where h{(j)) is the constant of integration. 

=^ c{9,ct>) = /i(</.)[sec^]3xi-024xio--r= 

2 

Now for r = rwMAP = ro (^^^7^) ' = 7.139 x lO^^. Finally 

=^ c{e,^) = /i((^)[sec^]^-6«>^i°"' 

Using the Gn equation one gets 

c(0,0) = <?((/>) [sec0] 1-56x10"' 

Therefore 

c(0» = ,g(0)[sec0]i-5'5xio-^ ^ h{cj))[sece]^-^^'''°~' 

Now in obtaining c{9, 0) we have taken « and neglected this term. So this expression is valid 

for those values of 9 where ^ g/a''^'' ~ i.e ^'^g^^ w constant and any deviation will be linear in 
9. Inserting the solution we find that the first derivative is proportional approximately to tan{9) 
equals 9 for small 9. This happens for those values of 9 for which c{9, (j))tan{9) « constant, which 
is true for values of 9 close to 0. Hence the given function describes anisotropy the best for values 
of 9 around 0. 

To evaluate the functions of (j) one must need another equation in c{9,(p) v/s (f> for which 
one must take higher order perturbations. But as of now b{t,9,(p) = g{(p)[sec9]^-^^^^^ ^t~^ = 
/i(^)[sec^]^ 6^xio ^— rpj^g diagonal terms of Gtj were set to zero. But since we now have the 
expression for b(t, 9, (p) we can now substitute it in the off diagonal terms say G12 and G21 and 
since G12 = 87rGTi2 and G21 = &itGT2i we have now actually got two terms for our stress energy 
tensor, i.e 

i - - 2 

3 



-g{4i){asec°'9tan9)tlmap 
-112 = J-21 = 



twmap 



to 
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where a = 1.56 x 10 ^ . Similarly the other components of the stress energy tensor can be 
determined. They would have derivatives of h{(t)). T12 and T21 involve h{4i). But it can be 
determined by not neglecting Thus the stress energy tensor can be completely determined and 
seen to correspond to inhomogeneity. This can be interpreted as introducing perturbations in the 
ideal fluid contained in Universe to obtain the non-ideal case. 



20 case 2 

In the above calculation we have pi non zero for first order and p2 = 0. But p2 = Goo for second 
order. Here 



P = 



6.765 X 10'** X (3 X 10**) 3 c((9, ^!>) + Kt' 



where po = 2.11 x 10 ^ and K is a constant t = twmap and to = presenttime 



Goo 





r 


4^2 




{toj 




34 





icr t3 



2cot{6)t'^ 



to 



de 



= Pl+P2 = 



Potp 



6.765 X 10'** X (3 X 10**)3c(6l, f— + Kr 



which yields 



{com X lO^)^^^^ - 'A^^«o 1^l2^ ^ t^,-2 ^ , ^ in-20 
06 



c{0, (/))(69 X 10'^) + Kt-'' + 4.35 x 10" 







This equation can again be solved for c{9, (p) , and hence for b and the complete Einstein tensor 
and stress energy tensor can be determined. The model is solved at WMAP time. 



21 Conclusion 

The correlation between the 2nd order gravity wave analysis which gives the space time 'wrinkles' 
on the FRW background and the density and curvature perturbations of second order that give 
the matter distribution on FRW background is expected, as the Einstein equations remain valid 
order by order. The correlation of the observed WMAP radiation data to the theoretical model 
is also suggested. The WMAP data has been explained by other models such as inflation driven 
inhomogeneity. In this work the contribution of gravity waves , non linearly coupled on the FRW 
background, organizing the density and curvature inhomogeneity was investigated. 
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